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This work is devoted to Bloch oscillations (BO) of cold neutral atoms in optical lattices. 
After a general introduction to the phenomenon of BO and its realization in optical lattices, 
we study different extentions of this problem, which account for recent developments in this 
field. These are two-dimensional BO, decoherence of BO, and BO in correlated systems. 
Although these problems are discussed in relation to the system of cold atoms in optical 
lattices, many of the results are of general validity and can be well applied to other systems 
showing the phenomenon of BO. 



I. INTRODUCTION 

This review is mainly addressed to researchers, who's prime scientific interests are far from the 
topic announced in the title. We assume the reader to have no preliminary knowledge of Bloch 
oscillations and introduce the problem step by step, beginning from the notion of optical lattices. 
On the other side, we try to avoid any extended derivations and the theoretical analysis of the 
considered phenomena is presented in a simplified form. Moreover, in some cases we only explain 
the main idea of the analytical approach (referring to the original papers) and directly proceed 
with the results. In this sense, this review serves only as an introduction to Bloch dynamics of cold 
atoms. For those already familiar with the subject we advise to skip the first sections and move 
directly to Sec. II VI and Sec. \V\ where the most recent developments in the field are discussed. 

A. Brief historical review 

Originally the problem of Bloch oscillations (BO) was formulated in context of crystalline elec- 
trons. Considering the response of the system to a static electric field, Zener came to the conclusion 
that, instead of the uniform motion which one would naively expects, the electrons in the crystal 
should oscillate pj. The period of these oscillations, now known as the Bloch period, is given by 
Tb = 2irfi/dF, where d is the lattice period and F is the magnitude of the static electric force. 
However, for a realistic strength of the electric field, the period Tb appears to be much smaller than 
the characteristic relaxation time r in the system. (The main contributors to r are the scattering 
by impurities or phonons and the electron-electron interactions.) Because of this reason BO have 
never been observed in the bulk crystal. 

The status of BO as a pure theoretical problem changed after fabrication of semiconductor su- 
perlattices 0|. Here, due to the essentially larger period d and the lower density of the carriers, one 
can satisfy the condition Tg < r and in 1992 the first experimental observation of BO was reported 
for such systems 0. It should be stressed, however, that BO in semiconductor superlattices are 
still dominated by the relaxation process. This difficulty is overcome by optical lattices, where 
standing laser waves and cold neutral atoms play the role of the crystal lattices and the electrons, 
respectively. In the latter system, the relaxation processes can be suppressed to any^desired level, 
which has offered unique opportunities for experimental studies of BO 0, E3, El V 



'Electronic address: kolovsky@mpipks-dresden.mpg.de URL: www.mpipks-dresden.mpg.de/~kolovsky 
^Electronic address: korsch@physik.uni-kl.de 



2 



We would like also to note that the semiconductor and optical lattices are not the only systems 
showing BO. As will be shown below, the deep origin of BO lays in the band spectrum of the 
system. In this sense, any spatially periodic system may show BO. A recent example is the 
periodic oscillation of a light beam in periodic photonic structures [Tel. fTll] . 

B. List of notations 

For the sake of quick reference we list below the notations used throughout the paper. 

• A - the laser wave length, defines the period of the optical lattices d = A/2; 

• pi = 2hki - double recoil momentum (kz, = 2tt/X); 

• E R = h 2 k 2 L /2M - recoil energy defines the characteristic energy scale of the system; 

• 4>a,K,(x) - Bloch states, i.e. the eigenstates of the quantum particle in a periodic potential; 

• E a (K) - energy spectrum of Bloch waves, with a being the band index and k the quasimo- 
mentum (—n/d < k < ir/d); 

• V'aiC 3 ') ~ Wannier states, which provide an alternative basis in the Hilbert space of the 
system, the site index I labels the wells of the optical lattice; 

• 9 at i(x) - Wannier-Stark states, i.e. the eigenstates of the quantum particle in a periodic 
potential plus a homogeneous field. Rigorously speaking, the Wannier-Stark states are res- 
onance or metastable states. 

• £a,i = E a i — iT a /2 - the spectrum of the Wannier-Stark states, where h/T a defines the 
lifetime of the states. 

II. OPTICAL LATTICES 

An optical lattice is a practically perfect periodic potential for atoms, produced by the inter- 
ference of two or more laser beams. In this section, we explain the origin of optical lattices, their 
properties and some limitations. 

A. Optical potential 

The physical origin of the optical lattice is the so-called dipole force, which acts on the atoms 
in the laser field. Indeed, let us consider a two-level atom in a standing plane wave: 

#= ( E I ) + m ( I l) ~ ^-Kt)cos(^x) ( J I) . (1) 

In Eq. (J2), M is the atomic mass, E g and E e are the ground and excited electronic states of the 
atom, Rabi frequency of the dipole transition between these states, ul the frequency, and kt the 
wave vector of the standing wave. Substituting the wave function ^(x, t) = exp(— iu}Lt)ip e (x, t)\e} + 
ip g {x,t)\g) into the Schrodinger equation with the Hamiltonian (JTJ and using the rotating wave 
approximation, one obtains the following system of coupled equations, 

ih ^d^' * = h5ip e (x,t) + ^i(; e (x,t) -hncos(k L x)^ g (x,t) (2) 
ih 9 ^ 9 ^^ = ^— ip g (x,t) - hQcos(k L x)ip e (x,t) , (3) 
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where 5 = (E e — E g )/Ti — ujl is the detuning. Let us now assume that the atom is initially 
in its ground electronic state and that the detuning 5 is much larger than the Rabi frequency 
f2. (More precisely, 5 is the largest characteristic frequency of the system.) Then ip e (x,t) ~ 
(17/(5) cos(kLx)ip g (x,t) and we end up with the Schrodinger equation 

V(x) = V cos 2 (A: L x) , V = -hn 2 /5 , (5) 

which describes the motion of the atom along the standing wave. The potential ©, which has a 
spatial period d equal to one half of the laser wave length, d = A/2, is called the optical potential 
or, simply, the optical lattice. Conveniently, the depth of the optical lattice is measured in units 
of the recoil energy Er = Ti 2 k\/2M . For example, for sodium atoms in a laser field detuned by 
60 GHz from the Z?2 sodium line (resonant wave length 589 nm), a 4 mW power laser creates an 
optical potential with an amplitude Vo of about 10 recoil energies |l2j]. 

B. Spontaneous emission 

In a more sophisticated approach, where the electromagnetic field is treated quantum- 
mechanically, the dipole force appears due to the stimulated exchange of photons between the modes 
of the electromagnetic field, associated with two counter-propagating running waves. Namely, the 
atom absorbs a photon from one of the running waves and 'immediately' emits it into the other 
wave, getting a recoil kick pi along the standing wave. During this absorption-emission process, 
the atom may emit a photon in the other modes of the electromagnetic field, getting a recoil kick 
in an arbitrary direction. The latter process, known as spontaneous emission (which should be 
opposed to the stimulated emission discussed above) is a kind of relaxation process due to interac- 
tion of the system with the environment, i.e. a bath of the electromagnetic modes. The rate 7 of 
spontaneous emission is given by the product of the natural width of the excited level 7 (which is 
a unique characteristic of the chosen atomic transition) and the population of the upper state, 

7 = ^(Q/5) 2 . (6) 

This equation implies that by simultaneously increasing the detuning and the intensity of the laser 
field, one can keep the depth of the optical potential constant but suppresses the interaction of the 
system with its environment. For example, in the case of sodium atoms in a laser field detuned by 
60 GHz, the rate of spontaneous emission is 7 ~ 100 s -1 , which is actually negligible on the time 
scale of the laboratory experiment [l2T |. 

C. Lattice dimensionality 

Up to now, we have considered an idealized situation of plane waves. In practice, however, one 
deals with beams of a finite width, i.e., 

cos 2 {k L x) , (7) 

where ro ~ 50 [im is the 1/e diameter of the beam. Note that for a red detuning, the laser field 
also provides a transverse confinement for the atoms. Besides the optical potential Q), there is 
an additional harmonic confinement Vt rap (r) ~ co 2 x 2 + uj 2 y 2 + oj 2 z 2 in the laboratory experiment 
due to a magnetic time-orbiting potential, which is used to capture the atoms during the cooling 
procedure. After the sample preparation (the sample preparation includes the cooling of the 
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FIG. 1: Energy spectrum of an atom in an optical potential of depth Vq = Er (left panel) and Vq — 4Er 
(right panel). 

atoms and an adiabatically switching on of the optical potential) this harmonic potential can be 
kept 'switched on' or relaxed towards zero. If not stated otherwise, we assume the second case 
throughout the paper. 

Using two crossed standing laser waves (4 running waves) one can create two-dimensional lattices 
with approximately (tq/X) ~ 100 wells in each direction. There is a high degree of freedom in 
choosing a particular form of the 2D potential. For example, by playing with the frequencies of the 
waves, one can realize separable 'egg-crate' or non-separable 'quantum dot' potentials (see figure |S] 
below); changing the angle between the crossing standing waves from 90 to 60 degrees transforms 
the square lattice into a hexagonal one; and so on. We shall study 2D optical lattices in more 
detail in Sec. II VI Needless to mention that using three mutually perpendicular standing waves one 
gets a true 3D lattice - an extention of the results of Sec. II VI to the 3D case is straightforward. 

Going ahead, we note that, beside providing a richer dynamics of BO, the higher dimensionality 
of an optical lattice also affects the strength of the atom-atom interaction. Namely, due to a stronger 
confinement of the atomic wave function, the effective constant of atom-atom interactions in 2D 
lattices is at least one order of magnitude (two - four orders for 3D lattices) larger than in ID 
lattices (see Sec. IV B|) , 



As well known, the energy spectrum of a quantum particle in a periodic potential consists of the 
Bloch bands. An example of the atomic Bloch band spectrum is given in Fig.^ where zero energy 
corresponds to the bottom of the potential wells. If the optical lattice is switched on adiabatically, 
the atoms populate the bottom of the ground band. The characteristic width of the distribution 
over the quasimomentum k depends mainly on the frequency of the harmonic trap. For a small 
frequency (a weak confinement) the atoms may coherently populate hundreds of the wells, which 
results in a very narrow distribution in the quasimomentum, Ak ~ O.Olir/d. Thus one may speak 



D. Bloch waves 
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of an atomic Bloch wave, 

<t>a,K,{v) = exp(ira)xa, K (x) , Xa, K ( x + d) = Xa,n(x) . (8) 

One of the facilities provided by optical lattices is that the atomic Bloch waves can be directly 
measured in the laboratory experiment ^3). Usually, the measurement goes as follows. After 
preparation of the Bloch wave, the optical potential is abruptly switched off and the atoms move 
in free space for a given time (the so-called 'time-of-flight'). Then the atoms are exposed to 
resonant light and, by 'taking a picture' of the atomic cloud, one records the spatial distribution 
of the atoms. Because the time-of-flight is known, this spatial distribution carries information 
about the momentum distribution of the atoms in the optical lattice. The latter is given by the 
squared Fourier transform of the Bloch wave @ and consists of a number of peaks, separated 
by PL = 2TT%/d. This peak-like structure of the momentum distribution (see Fig. [2] below), well 
observed in the laboratory experiments, is a direct indication of the atomic Bloch waves. 

We conclude this section by introducing the Wannier states, which we shall use later on. These 
Wannier states are obtained by integrating the Bloch states (JSJ) over the quasimomentum, 

/■n/d 
dKexp(-idKl)(j) aK (x) , (9) 
-TT/d 

and provide an alternative basis in the Hilbert space. Unlike the Bloch waves, the Wannier states 
are localized in space. Note that, because of = il> a ,o{x — Id), it suffices to calculate only one 

Wannier state (/ = in what follows) for each energy band. 



III. BLOCH OSCILLATIONS IN ID LATTICES 



This section studies different regimes of BO of cold atoms in ID optical lattices. It is implic- 
itly assumed in what follows that neither spontaneous emission nor atom-atom interaction are 
important and, thus, we can use the single-particle Schrodinger equation to analyze the problem. 
Beside this, we assume that the transverse motion of the atoms is frozen (i.e., we are dealing with 
a ID problem). Although we do not discuss the validity of this approximation, the experimental 
results jj, [5|, indicate that this is, indeed, the case realized in quasi ID lattices. 



A. Bloch period and Landau-Zener tunneling 



Bloch oscillations are the dynamical response of the system to a static force: 

-2 

H = H + Fx, H = ^ + V Q cos 2 {k L x) . (10) 

For neutral atoms, the static force F is usually introduced by accelerating the optical lattice, 
V(x) -» V(x - at 2 /2), which can be done by an appropriate chirping of the frequencies of two 
counter-propagating waves. Then, in the lattice coordinate frame, the atoms experience an inertial 
force of magnitude F = aM. The other option is to employ the gravitational force for a vertically 
oriented optical lattice, then a = 9.8 m/s 2 . 

The common approach to the problem of BO is to look for the solution as a superposition of 
Houston functions [13], 

^(x, *) = c a(t)^a(x, t) , (11) 
a 



ip a (x,t) = exp (~^J dt'E a (K')j <t> a , K t{x) , 



(12) 
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FIG. 2: Dynamics of the atomic momentum distribution P(p), induced by a weak static force F < 0. The 
depth of the optical potential is Vb = IOEr. The weak wiggling of P(p) is an artifact due to the finite size 
of the lattice (L = 10, periodic boundary conditions) used in the numerical simulations. 

where <t> a n'{ x ) 1S the Bloch function with quasimomentum k' evolving according to the classical 
equation of motion p = —F, i.e k' = kq — Ft/%. Substituting Eq. (fT2|) into the time-dependent 
Schrodinger equation with the Hamiltonian (jlUj) . we obtain 

ihc a = F^X Qi/3 ( K ')exp (~ J* dt'[E a {K') -^(k')]) cp , (13) 

where X ai p(K) = J dxXa K (x)dxf3, K ( x )/dK. When neglecting the inter-band coupling, i.e X a ^{n) = 
for a 7^ (3, we have i%c a = FX a ^ a (K')c a and, thus, 

\Ca(t)\ = \c a (0)\ . (14) 

This solution is the essence of the so-called single-band approximation. 

The correction to the solution l|14|l is obtained by using the formalism of Landau-Zener tun- 
neling. In fact, when the quasimomentum k' explores the Brillouin zone, an adiabatic transition 
occurs at the points of 'avoided' crossings between the adjacent Bloch bands (see, for example, 
the avoided crossing between the 1st and 2nd band in Fig. ^a) at k = ir/d). As a result, the 
population of the ath band decreases exponentially with the decay time 

r = h/T a , T a = aaFexpi-ba/F) , (15) 

where a a and b a are band-dependent constants. Note, that Eq. ()15|) provides only an estimate 
for the mean decay rate and, to find the exact dependence T a (F), one has to employ a different 
approach which we shall briefly discuss later on in Sec. IIIIDI 

As follows from the estimate IJ15JI for a weak static force the Landau-Zener tunneling can be ne- 
glected and the solution of the problem is essentially given by Eq. (|12[). i.e., ip(x, t) ~ 4>a,K -Ft/h{ x )- 
The linear change of the quasimomentum of the Bloch wave results in a periodic change of the 
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atomic momentum distribution P{p) = \tp(p,t)\ 2 (see Fig. |2J) which is the quantity measured in in 
the laboratory experiments 0, H|. (Since the measurement is destructive, one has to repeat the 
experiment several times to record the time-evolution of the momentum distribution.) The period 
of these oscillations is given by the Bloch period, 

T B = 2-Kh/dF (16) 

with d = A/2 and F = aM. Using the single-band approximation, it is also easy to show that the 
mean momentum evolves as (p(t)) = Mv(Ft), where v(k) = dE a (n) /hdn is the group velocity. 
Note, that the amplitude of oscillations of (p(t)} is proportional to the band width and, thus, can 
be extremely small (deep optical lattices). The oscillations of the atomic momentum distribution, 
however, are qualitatively the same independent of the particular choice of the parameters and, in 
this sense, are a more reliable signature of BO. 

B. Wannier-Stark ladder 

Bloch oscillations can also be described in terms of the Wannier-Stark states, which provide a 
useful insight into the physics of this phenomenon. To simplify the analysis, we consider a tight- 
binding model - an additional (after the single-band) approximation to the original problem. This 
approximation is not crucial and one obtains similar results for the single-band model. Using the 
notion of the Wannier states Q, the Hamiltonian of the tight-binding model has the form, 

Htb = J2(E a + dFl) + + l )W + \ l ~ 1>M ) > ( 17 ) 

i z i 



where J a is the hopping matrix element, E a = E a (n), and (x\l) = tp a i(x). The Hamiltonian (|17jl 
can be easily diagonalized, giving the spectrum, 

E u (k) = E a + J a cos(cfo) , if F = (18) 

(note, in passing, that the tight-binding model approximates the Bloch dispersion relation by a 
cosine function), and 

E a>l = E a + dFl, if (19) 

The discrete spectrum P'Jfl is known as the Wannier-Stark ladder and the eigenfunctions corre- 
sponding to E a> i, 

|^> = £^(^)m, (20) 

(here J n {z) is the ordinary Bessel function) are known as the Wannier-Stark states. Because 
the Bessel functions J n { z ) are exponentially small for \n\ > \z\, the Wannier-Stark states are 
localized in space with a localization length ly/s = 1 ( m units of lattice period) for dF > J a and 
1 W S ~ J a /dF for dF < J a . One may refer to these two cases as the 'strong force' and 'weak force' 
regimes. In this paper, however, we reserve the term 'strong force' to static force magnitudes which 
break the single-band approximation (strong Landau-Zener tunneling). 

C. Wave packet dynamics 



In Sec. 1111 Al we have considered the case of a Bloch wave as an initial condition. It is also 
interesting to study a situation, where only a few wells of the optical potential are populated [ill ]. 
Then, along with the oscillations in momentum space, the atoms also oscillate in configuration 



8 




FIG. 3: Spatial oscillations of the localized wave-packet. The amplitude of the oscillations is given by the 
localization length of the Wannier-Stark states Iws ~ 1/F. 

space. This is illustrated in Fig. |3] which shows the dynamics of P(x) = \tp(x,t)\ 2 for a 'minimum 
uncertainty' wave packet. The amplitude of the oscillations is given by the localization length of 
the Wannier-Stark states. Indeed, since the general solution of the Schrodinger equation can be 
written as a sum over the Wannier-Stark states, 

ij){x, t)=J2 c i exp(-iE7«,it/ft) *cmO) . ( 21 ) 
l 

the localization length Isw defines the maximum distance where the wave-packet can move to. 

For the purpose of future use, we also display the solution (|21|) in the momentum representation. 
Expanding the initial state in terms of Bloch states, ip(x,0) = J dKg a (%n)(j) a>K {x), we introduce 
the envelope functions g a (p) = 9a{fi>K>)- Then, using the extended Brillouin zone representation, 
the solution ijj(p,t) can be represented as 

oo 

^(p,t) = ^exp(-iE a t/h) 9a(p + Ft + n PL )^ afi (p) (22) 

a n=—oo 

(in comparison with Eq. Q21|) here we also included the sum of a). Equation (|22j) is well suited for 
a numerical simulation of Bloch dynamics and has actually been used in Sec. IHI^l to illustrate the 
oscillations of the momentum distribution. 

D. Bloch oscillations for strong static forcing 

Above we have considered only the case of a weak static force, where the Landau-Zener tunneling 
is negligible and one can use the single-band approximation to study the system dynamics. This 
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FIG. 4: Imaginary part of the complex energies Il2.'ill as the function of the inverse scaled static force 1/F 
(F -> dF /2wV ). The depth of the optical potential is Vq = 8E R . 



section is devoted to the regime of strong forcing, which we shall analyze by using the formalism of 
metastable (or resonance) Wannier-Stark states. The latter are denned as non-hermitian eigenstates 
of the Hamiltonian (|10|). corresponding to the complex energies 

S a>l = E a + dFl - iT Q /2 . (23) 

The spectrum (j23|) generalizes the notion of the Wannier-Stark ladder ()19|) and is given by the 
poles of the scattering matrix of the system. For the details of the scattering matrix approach to 
the Wannier-Stark problem we refer the reader to the review Here we only note that this 

approach does not involve any approximation and, hence, the results presented below are rigorous. 

The decay constant T a in Eq. (|23|). which defines the lifetime of the Wannier-Stark states, has 
a rather nontrivial dependence on the static force. An example is given in Fig. |IJ where the six 
different curves correspond to Y a of the Wannier-Stark ladder, originating from the six lowest bands 
of the Hamiltonian Hq. Strong fluctuations of the decay rate, superimposed on the Landau-Zener 
dependence (|15|) . are noticed. These fluctuations are due to resonance tunneling, occurring when 
the positions of the Wannier-Stark levels in different wells of the optical potential coincide, i.e., 
when Re[£ a> i(F)] = Re[£^ m (F)}. 

Having the scattering problem solved, one generalizes Eq. (I22H for the system dynamics simply 
by substituting the stationary Wannier-Stark states Q2L)[) by the metastable Wannier-Stark states, 
the real energy E a by the complex energy E a = E a — iT a /2, and multiplying the whole expression 
by the Heaviside step-function Q(p + Ft) which truncates the momentum distribution at p < —Ft 

m 

oo 

ij}{p, t) = 6(p + Ft) exp(-iE a t/h) £ g a (p + Ft + np L ) * a , (p) . (24) 

a n=—oo 

The characteristic feature of the resonance Wanier-Stark states ^^,0 is an exponentially growing 
tail for negative momentum, ^f a> i(p) ~ exp(T a p/hF). Then, as follows from Eq. (|24j). the solution 
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ijj(p,t) is essentially a sequence of wave packets, separated by a distance pi. In the coordinate 
representation, these equidistant sequence of the 'momentum' wave packets transforms into a train 
of 'coordinate' wave packets, distributed in space according to a square law. Thus, in the strong 
field regime, the atomic array acts as a matter laser, emitting one pulse of matter per Bloch period 

1 

E. Related problems 

In this subsection, we briefly discuss two important modifications of the problem of atomic 
BO. The first one deals with BO in the presence of harmonic confinement 0, i.e. V(x) = 
Vq cos 2 (kz J x) + Mw^ 2 /2. (The characteristic value of the frequency u x is a few Hz, which should 
be compared with the frequency of small atomic oscillations ~ ojr^Vo/Er) 1 / 2 of few kHz.) Let 
us consider the situation where the atoms are located far from the trap origin. In practice, this 
initial condition is realized by a sudden shifting of the center of the trap to a distance xq S> d 
||. Then the atoms locally feel a static force F = Mu 2 xq, which is one of the preconditions 
for BO. It should be noted, however, that the analogy with BO should be drawn here with some 
precautions. Indeed, using the tight-binding approximation, the Hamiltonian of the atom in the 
combined potential reads as 

i i 

with v = Muj 2 d 2 . This Hamiltonian formally corresponds to the Hamiltonian of the quantum 
pendulum, 

Heff = -^ + JaCO S e, (26) 

and, hence, the characteristic frequency of atomic oscillations is given by the frequency of the 
pendulum. The latter is known to have a rather nontrivial dependence on pendulum energy and 
can be approximated by a linear law only in the asymptotic region of large energies ^3]. Thus, to 
mimic BO (in a sense that the frequency of oscillations is inversely proportional to the local static 
field) one should satisfy the condition that the pendulum (|2*fj)l is well above its separatrix. 

The second problem we would like to mention deals with BO in a resonant or near-resonant 
laser field. In this case, the external degree of freedom of the atom (the motion of the center of 
mass) cannot be decoupled from its internal (electronic) degree of freedom and, thus, we have to 
solve the system of partial differential equations (|2I3|) (with a static term added) exactly, without 
adiabatic elimination of the upper state [l8|. It was found, in particular, that BO of the atoms 
in a resonant field enforce a kind of Rabi oscillations, where up to 90 percent of the atoms may 
appear in the excited electronic state. We note, however, that an experimental realization of this 
interesting regime of BO requires a very narrow width of the optical transition h'j < En. 

IV. BLOCH OSCILLATIONS IN 2D LATTICES 

In this section, we study Bloch oscillations in 2D optical lattices. We restrict ourselves to square 
lattices, created by laser beams of equal intensities. If the frequencies of the crossing standing waves 
also coincide, the optical potential is given by 

V(x, y) = V [cos{k L x) + cos(/c L y)] 2 , (27) 

as can be easily shown by repeating the derivation of Sec. 2.1. Note that the potential (|27|) is not 
separable. If, however, the frequencies of the waves are slightly mismatched, we obtain a separable 
potential 



V(x,y) = V [ cos 2 (k L x) + cos 2 ] {k L y)] 



(28) 
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Obviously, the property of separability can be attributed only to 2D or 3D potentials. The conse- 
quences of this property for the dynamics of BO is one of the main questions we address in this 
section. 




FIG. 5: Contour plot of the non-separable optical potential Ij27(l (left panel), and the separable potential 
(|28|l (right panel). Note that the primitive translation vectors for the lattice (|27|l are rotated by 45 degrees 
with respect to the laboratory coordinate system. 



A. Band spectrum 

We begin with the analysis of the Bloch band spectrum for the specified 2D lattices. To simplify 
the equations, we shall use a scaling where the coordinate is measured in units of the lattice period 
d and the time in periods of the recoil frequency (t — ► u)Rt). This scaling (which also involves a 
rotation of the coordinate system) leads to the dimensionless Hamiltonian 



n 



d 2 



d 2 



dx 2 dy 2 



+ cos x + cos y — e cos x cos y 



(29) 



where the only independent parameter is the the scaled Planck's constant Ti ~ (Er/Vo) 1 / 2 and the 
constant e equal to zero or ±1 for a separable and non-separable potential, respectively. (Although 
e can take only the specified integer values, for theoretical purposes it might be useful to consider 
the whole interval < |e| < 1.) 

To characterize the two-dimensional dispersion relation E a (n x , K y ), we consider the cross-section 



of the spectrum along the lines n v = for k x < and k v 



for k x > 0. The result is depicted 



in Fig. El where the left and right panels refer to the case e = and e = 1, respectively. It is 
seen in the figure that a nonzero value of e strongly modifies the central part of the spectrum. (In 
particular, it removes the degeneracy between the bands along the lines K y = ±k x .) As concerning 
the ground Bloch band, the difference shows up in the coefficients of the Fourier expansion of the 
dispersion relation 



Eq(k x , K y ) = 22 Jm,n exp(i27T?7lK x ) exp(i27T?7.Kj / 



(30) 
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In the separable case, only the coefficients J m ^ n with n = or m = differ from zero while in the 
non-separable case all elements have non-zero values. (Also note the symmetry J± n ,±m = Jm,n-) 
At the same time, even the largest non-trivial coefficient Ji.i is only 1/20 of the coefficient Jo,i, 
which practically alone determines the width of the ground Bloch band. 



B. Fractional Wannier-Stark ladder 



Let us briefly discuss the energy spectrum of the 2D Wannier-Stark system, 

H = H + F x x + F y y , (31) 

where Hq is given in Eq. (|29|) . First we consider the special case when the vector F of the static 
force is parallel to one of the crystallographic axes of the lattice (the x-axis, to be certain). In 
this case, the Hamiltonian (|31|) possesses the 'ladder symmetry' along the field direction and a 
translational symmetry in the direction perpendicular to the field. Thus, the spectrum of the 
system consists of replica of the Bloch band £ a (K y ), shifted relative to each other by the Stark 
energy 2ttF (here 2tt stands for the lattice period). 

The above result can be extended to the case of arbitrary 'rational' directions of the field, 



where q, r are co-prime integers. In this case, one uses the transformation of the coordinates 19], 

, _ qx + ry , qy - rx 

X ( r 2 +g 2)l/2> V ( r 2 + g 2)l/2' 
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which introduces a new lattice with the period d! = 2ir(r 2 + q 2 ) 1 / 2 and matches the vector F to 
the primitive vector of this new lattice, 

~' 2 "'2 

H' = ^- + ^- + V(x',y') + Fx' . (34) 

It is also easy to see that the transformation ()33|) actually introduces s = r 2 + q 2 different 
(sub)lattices, whose Hamiltonians differs from (|34[) by an additive term (d'F/s)j, j = 1, . . . ,s — 1. 
Thus, for rational directions of F, the spectrum of the original Hamiltonian ()31|) is a fractional 
Wannier-Stark ladder along the direction of the field, constructed from the Bloch bands with y^s 
times reduced Brillouin zone in the direction perpendicular to the static force, i.e., 

£a,l( K ±) = ^( K -L) + ^2^^1/2 ' £a{K ± ) = E a (K ± ) - i^-^^ . (35) 

Note that for a separable potential the sub-bands £ a (K±), — s -1 / 2 < k± < s -1 / 2 , have zero width 
for any direction of the fields 9 = arctan(r/q), except 9 = 0, tt/2. (For the imaginary part of the 
dispersion relation one obviously has T a (F, 9) = T a (F cos 9) +T a (F sin 9), where T a (F) is obtained 
by solving the ID problem.) As e deviates from zero, the sub-bands gain a small but finite width. 
A more dramatic consequence of the non-separability, however, is the strong dependence of the 
decay rate T £ An±) on the quasimomentum k±, where the decay rate may vary by several orders of 
magnitude |20l |. 

C. Wave packet dynamics 

It is interesting to compare the wave packet dynamics for separable and non-separable poten- 
tials. Restricting ourselves to the ground Bloch band, the results of these studies |2ll . l22j | can be 
summarized as follows. 

For a weak static force (negligible Landau-Zener tunneling) and e = 0, the two-dimensional BO 
are given by a superposition of the one-dimensional BO. In other words, BO is a (quasi)periodic 
process with two periods defined by the projections of the static force to the crystallographic axes 
of the lattice, T x<y = 2irfi/ dF X)V . In coordinate space, this is reflected in Lissajous-like trajectories 
of a localized wave packet in the xy-plane. Note that for e = 0, the motion of the wave packet is 
generally non-dispersive. (Exclusions are 9 = 0, 7r/2, where BO along one axis are accompanied 
by a dispersive spreading of the wave packet along the other axis.) A nonzero |e| < 1 only 
slightly modifies this dynamics which is actually not surprising, because the dispersion relation 
(|3*0|) for the ground Bloch band can be well approximated by a 'separable' dispersive relation, 
E(K x ,K y ) Ri (Jifi/2) cos(2itk x ) + (Jq^/2) cos(2iTKy). In terms of the energy spectrum (|33)) this 
approximation amounts to neglecting the width of the bands E a (n±). 

The case of a strong static force is essentially more complicated and here the difference between 
the separable and non-separable potentials appears on the qualitative level. Indeed, in the strong 
force regime, the atoms may escape out of the potential wells through a sequence of Landau-Zener 
tunneling transitions to higher bands. The details of this process crucially depend on the particular 
structure of the upper bands which, as seen in Fig. H3 is quite different for e = and e ^ 0. Similar 
to the ID case, one can analyze the tunneling of the atoms by using the formalism of metastable 
(now 2D) Wannier-Stark states. This analysis leads to the prediction that in the separable case 
the atoms may escape out of the potential wells only along the x- and y-axis of the lattice, while 
for the non-separable potential additional escape channels appear |2lJ. This is illustrated in Fig. [71 
which shows 'snapshots' of two-dimensional BO in the strong force regime for e = (left panel) 
and e = 1 (right panel). The additional chanel along the (l,l)-crystallographic direction is clearly 
seen in the right panel. (Also notice a rich interference structure between the channels, which is 
absent in the separable case.) 
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FIG. 7: A fragment of the wave function of the atom in the separable (left panel) and non-separable (right 
panel) potential in the case of a strong static force. The system parameters are Ti = 2, F — 0.2, and 
FxlF y = 1. 

D. Related problems 

The 2D optical lattice also offers an opportunity for studying a number of related problems like, 
for example, 'two-bands' BO. Indeed, superimposing the potentials shown in Fig. one can easily 
realize the case, where the two lowest Bloch bands of the combined potential are separated by an 
arbitrary small gap (located along the edges of the Brillouin zone). Then, performing a BO, the 
atom tunnels between these two bands at each crossing of the Brillouin zone, which results in a 
very non-trivial dynamics of BO |2ll. l22l|. 

The other problem we would like to mention is the scattering of an atomic beam by a 2D 
optical potential [23j. In classical dynamics, and for a non-separable potential, this would be a 
chaotic scattering process, with fractal basins for the scattered channels. The classical and quantum 
scattering of the atoms by 2D lattices is studied in some details in the paper cited above. 

V. DECOHERENCE OF BLOCH OSCILLATIONS 

Since BO is a coherent quantum phenomenon, it is important to study the processes which 
cause a decoherence of the system. In this section we consider two of these processes - decoherence 
due to spontaneous emission and due to atom-atom interactions. We would like to stress that 
decoherence or relaxation, usually considered as 'unwelcome' phenomena, are also of interest on 
their own. Indeed, as was already noticed by Esaki [24|, the conventional conductivity is an 
interplay between BO and relaxation processes. Thus, the analysis of decoherence of BO is a 
necessary step for developing a theory of atomic conductivity. 

A. Decoherence by spontaneous emission 

As mentioned above in Sec. Ill Bl spontaneous emission is a particular case of the systems 
interaction with an environment. For this kind of problem, the approach based on the single- 
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FIG. 8: Dispersion (Ax 2 ) = (x 2 ) — (x) 2 of the atomic wave packet as a function of time. Time is measured 
in the units of Tj = 2Trh/ J, the scaled rate of spontaneous emission hj/ J — 0.05, d — it, and the scaled 
magnitude of the static force (F — ► dF/J) is indicated in the figure. The slopes of the dashed lines are given 
by the values of the diffusion coefficient i|39[l . 

particle Schrodinger equation is not applicable, and the dynamics of the system should be described 
in terms of the density matrix p(x,x';t), which is defined as the trace of the total wave function 
(system plus environment) over irrelevant variables of the environment. For a specified environment 
(the photon bath), the density matrix obeys the master equation [25| . 



where H is the Hamiltonian (|1U|) (we consider a quasi one-dimensional lattice), 7 the rate of 
spontaneous emission ©, and L u the projection of the recoil operator on the x axis, 



Note that Eq. (|36j) has the Lindblad form and, thus, Tr[p(i)] = J dxp(x, x;t) = 1. The distribution 
P(u) of the random variable u in Eq. (I36|) is denned by the angle distribution for the momentum 
of the spontaneously emitted photons |2a |. and is, in the case of linearly polarized light considered 
here, approximately given by P{u) = 1/2. 

Using the tight-binding approximation, i.e. substituting the Hamiltonian H by the tight- 
binding Hamiltonian (|17|) . and the recoil operator (|37|) in its tight-binding version L u = 
1)' exTp(i7rul)\l)(l\, the master equation (|36|) can be solved analytically, with the following 
main results 27]: The spontaneous emission leads to decoherence of the system, i.e. the density 
matrix tends to a diagonal one in the basis of the Wannier states. As a consequence, BO decay as 



dp 
dt 




(36) 



L u = cos(kLx) exp(iukLx) 



u\ < 1 . 



(37) 



(p(t)) = p exp(-z/t) sm(u) B t) 



(38) 



where the decay rate v appears to coincide with the rate 7 of the spontaneous emission. The decay 
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of BO is accompanied by (asymptotically) diffusive spreading of the atoms (see Fig. |HJ) 



{x\t))~Dt, S=(g) 




(39) 



(here po is the amplitude of BO in the absence of the relaxation process). Note that (since ujb ~ 
F) the static force actually suppresses the diffusion. Considering the diffusion as a 'generalized 
conductivity', this result agrees with the prediction of Esaki and Tsu that the conductivity of 
the system tends to zero when the frequency of BO becomes much larger than the characteristic 
frequency of the relaxation processes. 

The (numerical) analysis of the system dynamics beyond the tight-binding approximation leads 
to qualitatively the same results, although the quantitative deviation can be larger than 50 percent 

. It should also be mentioned that the validity of Eq. i(3*H|) still assumes a low population of the 
excited electronic state of the atom and, hence, the case of resonant driving (briefly analyzed in 
Sec. IIII El for 7 = 0) is excluded from this consideration. 



Up to now, we have studied BO using single-particle quantum mechanics. This is justified only 
for a very dilute gas of atoms, where the atom-atom interactions can be neglected. If this is not 
the case, the problem of BO becomes very diverse and, first of all, one should distinguish between 
Bose and Fermi statistics. In this review we restrict ourselves to bosonic atoms. Moreover, we 
assume in what follows that the initial state of the system is a Bose-Einstein condensate, where 
all atoms occupy the zero quasimomentum state of the ground Bloch band. To study the time 
evolution of this state, one usually uses the Gross-Pitaevskii (or nonlinear Schrodinger) equation. 
It is understood, however, that this equation has a limited applicability - in particular, the Gross- 
Pitaevskii equation is unable to describe the decoherence of the system. Because of this, we employ 
here a more general approach, based on the Bose-Hubbard model. 

The Bose-Hubbard model (with a static term added) generalizes the tight-binding Hamiltonian 
(|17|) to the multi-particle case. To simplify the analysis we shall consider only the ID Bose-Hubbard 
model, even when discussing the 3D lattices. (This approximation is not crucial for the phenomena 
discussed below.) Then the Hamiltonian of the system has the form 



(the band index a = is omitted, J > 0). In the Hamiltonian (|40j) . the creation operator aj and 
the annihilation operator ai 'creates' or 'annihilates' an atom in the Zth well of the optical potential 
in the Wannier state ipi{x). Hence, the first term on the right hand side of Eq. 1)40(1 is responsible 
for the tunneling (hopping) of the atoms between the wells of the optical potential. Note that the 
hopping matrix elements J (defined by the overlap integral of the Wannier states in neighboring 
wells) exponentially depends on the depth of the optical potential and, hence, can be easily varied 
by several order of magnitude (0.30 < J /Er < 0.0035 for 2 < Vq/Er < 22). The second term in 
the Hamiltonian (|40|) is the Stark energy of the atoms in the homogeneous field. The third term is 
the interaction energy of the atoms sharing one and the same well, where the interaction constant 
W is mainly defined by the s-wave scattering length a sc of the atoms and by the geometry of the 
lattice, 



As an estimate, one can use W = 0.28-Er - the experimental value for 87 Rb atoms (a sc = 5.8 nm) 
in 3D separable potential of the depth Vq = 22Er 28]. Interpolating this result to Vq = 2Er would 
give W = 0.027£j?. 



B. Interacting atoms and the Bose-Hubbard model 




(40) 




(41) 
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A remark about the phase diagram of the Bose-Hubbard model (the Hamiltonian (|4U[) without 
the static term) is appropriate here. As it is known, the Bose-Hubbard model shows the super- 
fluid/Mott-insulator quantum phase transition when the ratio between the parameters J and W 
exceeds some critical value 29]. (In the laboratory experiments one varies the ratio W/J by 
changing the depth of the optical potential - the critical value for 87 Rb atoms is Vq ~ 15-Er and 
Vq ~ 44.Er for the 3D and ID lattices, respectively [3CJ, |3J.) Obviously, the dynamical response 
of the system to a static force also depends on whether the system is in the super-fluid or Mott- 
insulator regime. In the next section, we study the Bloch dynamics in the super-fluid regime, where 
the ground state of the Bose-Hubbard system (which we take as the initial state in our numerical 
simulations) can be well approximated by the product of Bloch waves with zero quasimomentum. 
The response of the system to a static force in the Mott-insulator regime will be briefly discussed 
in Sec. IVDl 



C. Decoherence due to the atom-atom interactions 



Numerically, the problem of BO of interacting bosonic atoms consists of solving the 
Schrodinger equation for the multi-particle wave function ^f(t) = Xm c n(£)l n ); where |n) = 
| . . . , rfy_i, ni, ni+i, . . .) are the Fock (number) states, given by the symmetrized product of Wannier 
states ipi(x). Note that the translation symmetry of the system, obviously broken by the static 
term, can be actually recovered by using the gauge transformation, 

H - H(t) = -1 (j^e-^al^ + h.cj + y 5>(»*J - 1) . (42) 

This allows us to impose periodic boundary conditions, which greatly facilitate the convergence in 
the thermodynamic limit N, L —* oo, N/L = n. Here L is the lattice size and N the number of 
atoms. For given L and N, the dimension of the Hilbert space (the total number of Fock states) 




FIG. 9: Decay of BO due to the interaction induced decoherence. Parameters are J = 0.038-Er, W = 
0.032E R , and dF = -0ME R . 
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FIG. 10: Instantaneous spectrum of the Hamiltonian (|42|l . Lattice size L = 5, number of the atoms N = 5, 
interaction constant W = 0.1 J. 

is M = (N + L — 1)\/N\{L — 1)!. Note that the dimension of the Hilbert space but not the size of 
the system controls the convergence in the thermodynamic limit. 

Having a (numerical) solution of the Schrodinger equation, ihd\^ (t)) / dt = H(t)\^(t)), we then 
calculate the single-particle density matrix, 

p(x,x';t) = J2Mx)Mx')pi,m(t) > Pl,m{t) = {^(t)\aja m \^(t)) , (43) 

£,771 

which carries essential (although not complete) information about the system. In particular, the 
diagonal elements of the density matrix (|43j) in the momentum representation define the momentum 
distribution of the atoms P(p) which, as discussed above in Sec. Ill Dl is the quantity most easily 
measured in the laboratory experiments. 

An example of the time evolution of the momentum distribution is depicted in Fig.^JJ Comparing 
this figure with Fig.[2]for BO of non-interacting atoms, a rapid decay of BO is noticed. The reason 
for this decay is the decoherence of the density matrix, similar to that considered in Sec. IV Al but 
with the fundamental difference that here the system itself plays the role of a 'bath'. To get a 
qualitative understanding of this phenomenon, it is useful to consider the instantaneous spectrum 
of the time-dependent Hamiltonian ()42|) (see Fig. HUj) . The thin line in Fig. ^TOJ is the mean-field 
solution (diabatic continuation of the ground state), where all atoms would oscillate coherently. 
However, due to Landau-Zener transitions at avoided crossings, the other Fock states become 
populated when the static force drives the system along the the mean-field solution. This leads to 
a thermalization of the atoms and, as a consequence, to the decay of BO. 

In principle, one can try to describe the thermalization process by thoroughly analyzing the 
Landau-Zener tunneling at the avoided crossings. An alternative (and, actually, more constructive) 
approach is to study the properties of the Floquet-Bloch operator 32], which we define as the 
evolution operator over one Bloch period: 
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FIG. 11: Quasiperiodic Bloch oscillations. Parameters J and W are the same as in Fig. but dF = Er. 
Note that the scale of the time axis is changed in comparison with Fig. 

It has been found that for the parameter region of a typical experiment with cold atoms in the 3D 
lattices (W ~ J, n ~ 1, and dF < J) the matrix of the Floquet-Bloch operator (|44|) can be well 
identified with a random matrix of the circular orthogonal ensemble and, thus, the system 1)40(1 is 
a quantum chaotic system [331. l34|. In fact, this is precisely the Quantum Chaos, which justifies 
the use of the terms 'bath' and 'thermalization' for the system of interacting atoms. 

D. Different regimes of BO of interacting atoms 

As already mentioned in Sec. IV Bl BO of interacting atoms is a rather diverse problem and the 
irreversible decay of BO discussed above is in no way the only possible regime of Bloch dynamics. In 
particular, a static force with magnitude dF ^> J suppresses the Quantum Chaos and BO become 
quasiperiodic (see Fig. Ill|). This regime of BO can be treated analytically and, for example, for 
the mean atomic momentum we have [35| . 

{p(t)} = po sm(uBt) exp{— 2n[l — cos(cjjy^)]} > = W/Ti . (45) 

Note, in passing, that the same interaction frequency uw = W/H (but for a different problem) has 
recently been observed in a laboratory experiment [2a |. 

The other limiting case of BO is the case of a weak atom-atom interactions, which is typically 
realized in the quasi ID lattices. Indeed, because of a weak transverse confinement ro 3> A, the 
integral J ipf(r)dr in Eq. (|41j) is orders of magnitude smaller than for 3D lattices. An increase of 
the filling factor from n ~ 1 in the 3D case to n ~ 100 in the ID case (present days situation) 
may not be able to compensate this decrease of the interaction constant W and, therefore, a 
regular dynamics of BO can be expected. Let us also remind that a large filling factor and a 
small interaction constant are usually considered as a validity condition for the Gross-Pitaevskii 
equation. We shall discuss this issue in more detail in Sec. IV1I 

It is also interesting to study 'Bloch oscillations' in the Mott-insulator regime ( J <C W). In this 
case, a response of the system to a static force has a resonant character (30, with the main 
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resonance corresponding to the condition dF ~ W. Providing resonant (or near resonant) forcing, 
the particle-hole excitations of the Mott-insulator states are created dynamically. This excitation 
process is reflected in the (almost) periodic dynamics of the atomic momentum distribution 37]. 
It should be stressed, however, that in spite of the formal analogy with BO, the origin of these 
oscillations as well as the characteristic period Tj ~ fij J are fundamentally different. 

VI. CONCLUSION 

We considered the phenomenon of BO for cold neutral atoms in optical lattices. In Sec. Inland 
Sec. Mil of the review, which mainly serve as tutorials, we focused on a dilute gas of atoms in 
quasi ID lattices. In addition to the theoretical analysis, we also discussed the scheme of a typical 
laboratory experiment on BO and indicated the characteristic values of the parameters. 

It can be safely stated that in the above case of a dilute gas in ID lattices the BO of cold atoms 
are well understood. The further progress in the field is related to the problems of BO in lattices 
of higher dimensionality, BO in the presence of relaxation processes, and BO of interacting atoms. 

A dilute gas of cold atoms in 2D lattices was discussed in Sec. IIVI Naively, one may expect 
the two-dimensional BO to be a superposition of one-dimensional BO, where the values of the 
static force F Xj y are given by the projection of the vector F to the crystallographic axis of the 
lattice. However, this is true only in the separable case. For a non-separable optical potential, the 
two-dimensional BO generally cannot be expressed in terms of the ID problem. This is especially 
the case for a strong static force (strong Landau-Zener tunneling), where non-separability of the 
potential manifests itself in a number of effects. 

The ultimate goal for studying BO in the presence of relaxation processes is to obtain a directed 
diffusive current of atoms, similar to that of electrons in a conductor. In Sec. II V A1 we studied the 
relaxation (decoherence) of BO due to spontaneous emission. Note, that the rate of spontaneous 
emission can be increased to any desired level by simply tuning the laser frequency closer to 
the atomic resonance. The effect of spontaneous emission is shown to lead only to 'non-directed' 
diffusion. The theoretical analysis of the other relaxation mechanisms, like scattering on 'impurities' 
is strongly in need. In this connection we would like to mention the recent experiment which 
studies the dynamics of the fermionic 40 K atoms with an admixture of bosonic 87 Rb atoms in 
inhomogeneous optical lattices. 

Finally, we comment on the problem of BO for a Bose-Einstein condensate. This field of 
research opens unlimited perspectives for studying the different phenomena of correlated systems. 
In particular, by loading a BEC of cold atoms into a quasi ID lattice, one can realize the case of a 
large filling factor (mean number of atoms per lattice cite) and a small interaction constant. This 
is believed to be the realm of the Gross-Pitaevskii equation, where BO becomes a macroscopic 
quantum phenomenon. We intentionally did not discuss this case because even the mentioning of 
all publications on this subject would draw this review out of the length limit. Instead, we analyzed 
BO of interacting atoms in the 3D lattice. Loading BEC in the 3D lattice increases the interaction 
constant by several orders of magnitude and simultaneously decreases the filling factor to n ~ 1. 
In this regime the Gross-Pitaevskii equation fails to describe the dynamics of the atoms and a 
fully microscopic treatment of the system is required. Such a microscopic approach is provided 
by the Bose-Hubbard model. It was shown that, depending on the value of the static force, BO 
of interacting atoms may vary from quasiperiodic to irreversibly decaying behavior. Moreover, in 
the latter case, the system appears to be chaotic in the sense of Quantum Chaos. This intriguing 
regime of BO is waiting for experimental studies. 

Financial support by Deutsche Forschungsgemeinschaft (Grand No. 436 RUS 113/658/0-1) is 
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